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THE ERDOS-SZUSZ-TURAN DISTRIBUTION FOR 
EQUIVARIANT PROCESSES 

JAYADEV S. ATHREYA AND ANISH GHOSH 


Abstract. We resolve problems posed by Kesten and Erdos-Sziisz- 
Turan on probabilistic Diophantine approximation via methods of ho¬ 
mogeneous dynamics. Our methods allows us to generalize the problem 
to the setting of general measure-valued processes in R n , and obtain 
applications to the distribution of point sets which occur in higher¬ 
dimensional Diophantine approximation and the geometry of translation 
surfaces. 
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1. Introduction 


1.1. Dirichlet’s Theorem. The most classical result in Diophantine approxima¬ 
tion is Dirichlet’s Theorem, which is stated in two forms: First, for all irrational 
a £ [0,1], and any integer Q > 1, there is a 1 < q < Q and a p £ Z relatively prime 
to q so that 

(Dl) \aq-p\< 2. 


As a corollary, one obtains that there are infinitely many ^ £ Q (with gcd(p , q ) 
satisfying 


(D2) 
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In this paper, we consider the distribution of the number of solutions to modified 
versions of m and (1D2I) from a probabilistic perspective initiated by Erdos-Sziisz- 
Turan m and Kesten [53] . Our methods show how this problem can be generalized 
and solved in many other geometric and number theoretic settings. 


1.2. The Erdos-Sziisz-Turan and Kesten distributions. In 1958 m , Erdos- 
Sziisz-Turan introduced a problem in probabilistic Diophantine approximation: 
what is the probability f(N,A,c) that a point a chosen from the uniform dis¬ 
tribution on [0,1] has a solution | G Q to the modified Diric.hlet equation 


( 1 . 1 ) 



with denominator q G [TV, ciV]? Here A > 0,c > 1 are fixed positive parameters, 
and N is a parameter which goes to infinity. We note that by a well known result 
of Hurwitz, A = -^= is the best allowable constant so that HU has infinitely 
many solutions for all a. Let EST (A,c,N) denote the random variable counting 
the number of solutions to m- Then, the Erdos-Sziisz-Turan question is the 
existence of the limit 

lim P(EST(A, c, N) > 0). 

JV—>oo 

Considering analagously a modified version of m, Kesten [21 j defined the se¬ 
quence of random variables K(A, N ) (our notation differs from [24]) as the number 
of solutions to 

(1-2) \aq -p\ < 1 < q < N. 

We consider the following 


Question. What (if they exist) are the limiting distributions of the random vari¬ 
ables EST (A, c, N) andI\(A,N)? 


Our main results shows that the limiting distribution exists and can be viewed as 
the probability of a random unimodular lattice intersecting a certain fixed region. 
Let X 2 = SL(2, R)/SL(2,Z) denote the space of unimodular lattices in R 2 ,via the 
identification 

g SL(2, Z) i-a gl?. 

Let fi -2 denote the Haar probability measure on X 2 , and given A G X 2 , A = gl?, let 
Aprim be the set of primitive vectors in A. 

Theorem 1.1. The limiting distribution of the random valuables EST(A,c, N) and 
K(A,N) exist and denoting the random variables with these limiting distributions 
as EST(A, c) and K(A), we have 

(1.3) P(EST(A, c) = k) = (A G X 2 : #(A prim O H A ,c) = k), 
and 

(1.4) P(K(A) =k)= ^(A g X 2 : #(A prim n R A ) = k) 
where 

(1.5) H AjC = {(x,y) £S? : xy < A, 1 < y < c}, 
and 

( 1 . 6 ) 


R a = {(x, y) G R 2 : |x| < A, 0 < y < 1}. 
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Expressing the limiting distributions as distributions on the space of lattices al¬ 
lows us to apply classical results on the geometry of numbers to obtain moment 
and concentration estimates ('Theorem 12. 4[) . Furthermore, our translation of the 
problem into a geometric and dynamical problem is axiomatic and flexible, and as 
we will see in 21 applies to point sets associated to linear forms and translation 
surfaces. Our proof also shows that the same result holds as long as a is chosen 
from a probability measure with a continuous density (and in higher dimensions, 
certain natural classes of singular measures). Our main result (Theorem 14.11) de¬ 
scribes how to define, construct and compute related distributions in the general 
setting of equivariant processes, which we define precisely in (2) 

1.3. History. This circle of problems has a long history, starting with the original 
paper of Erdos-Sziisz-Turan m- There, in addition to posing the problem of the 
limiting probability of P(EST(5, c, N) > 0) (they denoted this putative limit as 
f(A,c)), they showed that for A < jzngz, that the limit existed and 

f(A,c) = ^Alogc. 

TT Z 

Subsequently Kesten [24] showed that the the limit exists under the assumption 
Ac < 1, and this assumption was removed by Kesten-Sos [25] - Explicit formulas 
for f(A,c ) were computed much more recently independently by Boca [10] and 
Xiong-Zalrarescu [55] , who also considered localizing a to smaller intervals. Our 
results offer a significantly more refined picture of EST(5, c) and K(A ), yielding 
the limiting distribution (not just the probability of positivity) of both random 
variables in both this setting and in a variety of other geometric and number theo¬ 
retic contexts, and our methods can localize a to even shrinking intervals fi ll. 5. 21) . 
Recently, related results have been studied from the point of view of deviation of 
ergodic averages for toral translations, see Dolgopyat-Fayad da ng for intrigu¬ 
ing distributional results from this perspective. Our perspective is inspired by the 
work of Marklof-Strombergsson EDE2], who gave several beautiful applications of 
homogeneous dynamics to understanding the fine scale statistics of point sets. 

1.4. Organization. In the remainder of this introduction, we prove Theorem ll.il 
using equidistribution of horocycles on the space X 2 . In 21 we describe our general 
results in the settings of lattices; linear forms; diopliantine approximation on curves, 
and translation surfaces. In 21 we describe our general philosophy and state our 
main axiomatic theorem. In 21 we prove our results in their various incarnations 
on the space of lattices, and in 21 we prove our results in the setting of translation 
surfaces. In 21 we remark on how our approach can be applied in even more sit¬ 
uations, including (but not limited to) multiplicative diophantine approximation; 
complex diophantine approximation; the distribution of cut-and-project quasicrys¬ 
tals; and the distribution of discrete lattice orbits on the Clifford plane and related 
problems on cusp excursions on hyperbolic manifolds. 

1.5. Equidistribution on the modular surface. We relate the Erdos-Sziisz- 
Turan and Kesten distributions to dynamics on the space of unimodular lattices, 
and prove Theorem 11.11 This approach will generalize to higher dimensions, and 
will allow us derive a number of results using appropriate equidistribution results. 
We begin with the original question of Erdos-Sziisz-Turan and Kesten, proving 
Theorem 11.11 
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1.5.1. Proof of Theorem Al.il We note EST(A, c, N) = k if and only if 


3 exactly k distinct — £ 

q 


such that N < q < cN and 


a — 


< 


A 


Equivalently, there are exactly k vectors 


where 


€ Zp rim such that 


1 —a 

0 1 


G 9log nP^A,ci 


and g t = 


0 e" 


This follows by rewriting 


P 

a - 

q 


< -r, N < q < cN 

q 2 


as 


and then as 


q\qa-p\ < A, 


\xy\ < A, N < y < cN. 

Thus, we are interested in the measure of the set of a € [0,1] satisfying 
(1.7) ^(^log iV^a^prim bl Ha,c) ~ k. 

Let Xk denote the indicator function of the set 

{Ael : # (A prim n h A iC ) = k}. 

To compute P(EST(A, c) = k), we are interested in the N —> oo behavior of 


( 1 . 8 ) 


/ Xk (ffiog iVa a Zip r j m )(ia:. 

J o 


Let t]n denote the measure da on the set {gi og Aru a Zp rim : 0 < a < 1}, so we can 
rewrite m as t]n{xk)- By Zagier’s equidistribution theorem {49j p. 279], we have 


Vn 


M2, 


as N —> oo, where the convergence is in the weak-* topology. The functions Xk 
can be approximated by continuous functions with compact support on the space 
of lattices Xi , so by a standard approximation argument we conclude 


P(EST(A, c, N) = k ) =t) N (xk ) —> M 2 (Xk)- 

N—> oo 

For the Kesten distribution, we note by a similar argument that K{A,N) = k if 
and only if 

ff (ffiog AW Q Zp r j m n Ra) = k. 

Thus, proceeding as above, we obtain 

□ 
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1.5.2. Measures and windows. The proof of Theorem II .11 in fact yields much more 
information. A strengthening of Zagier’s theorem due to Shah m allows us to 
obtain the equdistribution result for any absolutely continuous measure on [0,1]. 
Thus, the limiting random variables EST and K do not depend on the initial 
distribution of a (as long as it is continuous). A different strengthening of Zagier’s 
result is due to Hejhal 22] (with subsequent work of Strombergsson (42]), which 
implies that we can sample a from smaller subintervals depending on IV, as long 
as the subintervals shrink no faster than IV -1 / 2 . 

We can also consider the Erdos-Sziisz-Turan and Kesten distributions associated 
to solutions of ED and El with C\N < q < c 2 N , with 0 < Ci < c 2 . The limiting 
distributions will again be given by the probability random lattices intersect fixed 
regions in k points, with the regions being given by 

Ha, cue = {(x,y) e M 2 : xy < A, a < y < c 2 } 

and 

Ra, Ci ,c 2 = {{x,y) e R 2 : \x\ < A,d <y < C 2 }. 


1.5.3. Moments and concentration. To compute moments of the Erdos-Sziisz-Turan 
and Kesten distribution, we need to understand the quantities 

OO 

^fc‘P(X = k),t € M, 

k =0 

where X is either EST(A,c) or K(A). For t = 1, we can rewrite this as 


ix 2 


(Ap r i m n Ha, c ) dy 2(A). 


By the Siegel mean value theorem [40], we have 

[ # {A prim n H a , c } dy 2 (X) = \ \H A c \ = ^AlogC. 

Jx 2 7T 7T 

This is the expected number of solutions to El (in the N —> 00 limit). Note that 
if A < t^j, 

— 1 + C Z ’ 


12 


P(EST(A,c) > 0) = V P(EST(A,c) =k) = -^A\ogC , 

Z- J nrZ 


k>l 


so we have 

P(EST(A, c)) = y /fcP(EST(A, c) = k) = y P(EST(A, c) = k ), 
k >1 k >1 

so for k > 1, P(EST(A,c) = k) = 0, that is, the Erdos-Sziisz-Turan distribution 
is concentrated at 1. For the Kesten distribution, similar computations yield the 
mean, 

E(K(A)) = -\R a \ = — . 

7T Z 7T Z 

In the setting of higher-dimensional Diophantine approximation, we will obtain 
bounds on higher moments via classical results on the geometry of numbers. 
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2. Erdos-Szusz-Turan and Kesten distributions in higher dimensions 

2.1. Diophantine approximation. We start with a natural higher dimensional 
generalization of the original Erdos-Sziisz-Turan problem. Let d > 2 = m + 1 
(d = 2 ,m = 1 corresponds to our original problem), and fix A > 0,c > 1 and 
a norm || • || on R m . Let x be chosen from the uniform distribution on [0, l] m , 
and let EST<j(A, c, N ) denote the number of solutions (p, q) £ Z m x Z (with (p, q ) 
primitive) to the modified Dirichlet equation 

(2.1) IN - p|| < Aq-*, 

with q £ [N, cN], and Kd(A,N ) denote the number of solutions to 

(2.2) ||x<2 - pH < AN~i, 

with q £ [1, N], 

Theorem 2.1. The limiting distributions o/ESTd(A, c, N) and K d (A , N) exist and 
the distributions of the limiting random variables ESTd(A,c) and K d (A) are given 
by 

P(EST d (A, c)) = p d (A £ X d : # (A prim n H d , A , c ) = k) 

and 

P(Kd{A)) = Hd (A £ X d : # (A prim n l?d,A) = &) 

where 

(2.3) ifd,A,c = {(x, y) £ x K : ||x||y < A, 1 < y < c} 

and 

(2.4) = {(x, y) £ x I : ||x|| < A, 0 < y < 1} 

2.2. Linear Forms. Next, we consider systems of linear forms. Let d = m + n, 
m, n > 1, fix A > 0, c > 1, and norms || • || m and || • ||„ on K m and 1". We consider 
the set of m linear forms in n variables, parameterized by M mxn (JSl), the set of 
mxn real matrices. We identify M mx „(M) with K 1 ™". Let X be chosen from the 
uniform distribution on [0, l] mn . We define the random variable EST mxn (A, c, N) 
as the number of solutions (p, q) £ Z m x Z” (with (p, q) primitive) to the modified 
Dirichlet equation 

(2-5) ||Aq—p|| m < A||q||"-, 

with ||g|| n G [N, cN]. Similarly, we define K mxn (A, N) as the number of solutions 
to 

(2.6) II ATq — p|| m < A\N\~ ™ , 

with ||g|| n £ [: 1,N]. 

Theorem 2.2. The limiting distributions of the random variables EST mxra (A, c, N) 
and K mX n(A , N) exist and the distributions of the limiting random variables EST mx „(A, c) 
and K mxn (A) are given by 

P(EST mx „(A, c) =k) = nd{ A £ X d : # (A pr i m n Hmxn,A,c) = k) 

and 

P(K mxn (A) — /c) — hdi. A £ X d : ff (A pr j m D Rmxn,A ) = k) 


(2.7) Hmxn,A,c = {(x,y) £ x R” : ||x|| TO ||y|| n < A, 1 < ||y||„ < c}. 



THE ERDOS-SZUSZ-TURAN DISTRIBUTION 


7 


(2.8) Rmxn,A = {(x, y) e K m x r : ||x|| ro < A, 0 < ||y||„ < 1}. 

We note that Theorem [m is a special case of the above Theorem with n = 1. 

2.3. Approximation on Curves. The subject of metric Diophantine approxima¬ 
tion on manifolds studies typical Diophantine properties of points on manifolds. It 
is well known and easy to see using the Borel Cantelli Lemma, that almost every 
real number is not very well approximable. This means that the inequality 

\qx-p\ < l/\q\ 1+e 

has at most finitely many solutions. This result generalises easily to arbitrary di¬ 
mension. In 1932, K. Mahler conjectured that almost every point on the curve 
(a;, x 2 ,..., x n ) is not very well approximable. Mahler’s conjecture started the sub¬ 
ject and there have been many subsequent works, including recent dramatic ad¬ 
vances due to Kleinbock-Margulis, Beresnevich, Velani, and others. The constraint 
of lying on a manifold makes the subject considerably more complicated than clas¬ 
sical Diophantine approximation. Nevertheless, our approach can be used to com¬ 
pute Erdos-Sziisz-Turan and Kesten distributions for vectors lying on curves. Let 
d = n + 1 , and (f >: [a, b] —> R" be an analytic curve whose image is not contained in 
a proper affine subspace, and || • || denote a norm on R ra_1 . Let x be chosen from 
the uniform distribution on [a, b], and let EST^(A, c) denote the random variable 
counting solutions to 

|| q<fi(x) — p|| < Aq~d,N < q < cN. 

Let K^{A,N) denote the random variable counting solutions to 
|| q4>{x) — p|| < AN~^ , 1 < q < N. 

Theorem 2.3. The random variables EST^,(A, c, N) and A^(A,iV) have limit¬ 
ing distributions, and the limiting random variables EST^,(A, c) and K C / > (A) have 
distributions given by 

P(EST 0 (A,c) = fc)= p{AeX d : #(A prim n H d , A ,c) = k} , 

and 

P(/\ 0 (A, c) = k)= /i{Ae X d : #(A prim n R d ,A) = k} , 
where H d A,c and R d ,A are as in Theorem \2.1\ 

Remark: There is also an analogue of Theorem l2.2l for curves in the space of linear 
]R mrl , which is exactly parallel to Theorem 12.31 

2.4. Measures and windows. As in the setting of 1-dimensional approxima¬ 
tion, we can also work with q (or ||( 7 ||) in appropriate subranges of the form 
[c\N, C 2 IV] with appropriate changes to the limiting distributions (replacing the 
y range with [ci, C 2 ]). Additionally, choosing absolutely continuous measures also 
does not change the limiting distribution. 
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2.5. Moments and concentration. Classical results from the geometry of num¬ 
bers allow us to compute moments of the random variables EST and K. We recall 
the definition of the Siegel transform: given / € C^R") and A G Xd define 

/( A ) = /(*>)■ 

Aprim 

Siegel showed 

/ fdud = 7777 / f dm, 

Jx d C(«) Jr* 

where / is Lebesgue measure. Thus, the expectation of the random variables EST 
and I\ is given by times the volume of the regions HdA, c and Rd(A). Building 
on Siegel’s work, Rogers [35] and Schmidt m computed bounds for higher moments 
of /. These can be exploited to give precise moment estimates for EST and K , and 
yield non-trivial concentration phenomenon. For the variance, their results yield: 

Theorem 2.4. (Rogers [35} Theorem 4], Schmidt [HI Theorem 3], see also [8} 
Lemma 4.3],) Let X be either EST mxn (A,c) or K mxn (A), and let g = E(X). 
There is a constant Cd, depending only on dimension d = m + n, so that 

V(X) = E((X - gf) < C d g. 


In particular, for any T > 0, 

P{\X-g\>T s fg)< 


Proof. The second assertion is an immediate consequence of the first. For the first, 
we note both EST mxra and K m x n are random variables counting the number of 
lattice points in a bounded set ( H mxn {A,c ) and R mxn {A)). By [ 8 ] Lemma 4.3] 
(which is essentially contained in Rogers), we have that for any random variables 
of this type, 

E(X 2 ) <g 2 + C d g. 


In fact Cd can be chosen to be 8 f(d — 1 )/((d) (for d > 3). 


□ 


Remarks: It would be interesting to use Rogers’ [35] results on higher moments 
(bounds on L p norm for p < d) to obtain further concentration results. Note 
that since Rm.xn{A) grows polynomially in A and H mxn {A,c) polynomially in 
c and logarithmically in A, so do the expectations of EST and K , which gives 
explicit concentration bounds in terms of A and c. Finally, we note that Kesten [24[ 
Theorem 3] considered the d oo limit of Kd and proved Poisson behavior (under 
appropriate normalizations) using the method of moments. 

3. Translation Surfaces 

Our approach also yields information on the geometry of the set of holonomy 
vectors of saddle connections on translation surfaces. Given g > 1, an translation 
surface S of genus g is a pair S = (X, w), where X is a compact Riemann surface 
of genus g and w is a holomorphic 1 -form. A saddle connection 7 on S' is a geodesic 
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(in the flat metric determined by w) connecting two zeros of w, with none in its 
interior. The holonomy vector of 7 is defined by 

z 1 := / weC. 

•A 

The set 

As := {S 7 : 7 a saddle connection on S} 

is a discrete subset of R 2 with quadratic growth (cf. Masur (33]), that is there are 
constants 0 < Ci < C 2 so that 

C 1 R 2 < #(A S n 5(0, R)) < c 2 R 2 . 

We define the moduli space f l g of translation surfaces by considering equivalence 
classes of translation surfaces up to biholomorphism. This space is decomposed into 
strata R(a) consisting of holomorphic differentials with zeros of order aq,..., a*,, 
where a = (aq,..., ak) is an integer partition of 2 g — 2. Each stratum consists of 
at most 3 connected components [29], and there is a natural Lebesgue probability 
measure yu on the each stratum R(a), known as Masur-Veech measure. 

There is a natural SL(2, Reaction on the space which respects the decom¬ 
position into strata, and acts ergodically on each connected component of a stra¬ 
tum [341143] . The set As varies equivariantly under this action, that is 

AgS = 9^-S, 

where SL(2,R) acts on R 2 by the usual linear action. The fine-scale geometry of 
the sets As has been a subject of much recent investigation El [3lUS 07], and our 
approach allows us to define Erdos-Sziisz-Turan and Kesten distributions associated 
to translation surfaces. We note that for g = 1 , fii = X 2 , so this setting is another 
natural generalization of the original Erdos-Sziisz-Turan and Kesten problems. 

Let 9 £ [0, 27t) be chosen from the uniform distribution. We want to understand 
how well vectors in As approximate the direction 9 , in terms of their length. Given 
A > 0, c > 1, N > 0, define the random variables EST(S, N) and K(S , N) by 

EST(S,A0 =#(r 0 A s nH Ac ,N) 

and 

K(S, N) = # (rgAs fl Ra,n ), 

where 

Ha, c ,n = {A, y) e R 2 : xy < A, N < y < cN} 

and 

Ra,n = {(x, y) £ R 2 : \x\ < A,0 < y < N}. 

We say So £ fl g has circle limit measure y on f l g if the measures d9 on {gtrgS}o<e< 2 n 
converge to y. A result of Nevo [15] shows that for any stratum R, yu~ a.e. S' £'H 
has circle limit measure yn- 

Theorem 3.1. Suppose So £ has circle limit measure yo■ Then EST(So,A r ) 
and K(S 0 , N) both have limiting distributions, and denoting the random variables 
with this limiting distribution by EST(S) and K(S), we have 

P(EST(S 0 ) = k) = y 0 (S G fl g : # (A s n H A}C = k)) 

and 

P(K(So) = k) = yo (S £ n g : # (As fl Ra = k)). 
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In particular, for any stratum TL and pu~a.e. So £ H, 

P(EST(5'o) = k)=p n (SeH:# (A s n H A , C ) = k) 

and 

P(K(S 0 ) = k)= m (S£H:#(A s n R A ) = k). 

3.1. Lattice Surfaces. For particular highly symmetric surfaces, we can say more. 
We denote the stabilizer of the point S = ( X , oj) € f l g under the SL(2, R) action by 
SL(V, w). A translation surface S is called a lattice surface (also known as an Veech 
surface ) if SL(X, u>) is a lattice. The lattices that occur are always nonuniform, 
and the SL(2,R) orbit of S is closed, a copy of SL(2,R)/SL(V, w) in Ll g . For these 
surfaces, we have 

Theorem 3.2. Suppose So = (Xo, uiq) is a lattice surface, and write T = SL(Vo, wo). 
Let Mr denote the Haar probability measure on SL(2,R)/r. Then EST(So, V) and 
K(S 0 , N) both have limiting distributions, and denoting the random variables with 
this limiting distribution by EST(S'o) and K(So), we have 

P(EST(S 0 ) = k) = p r (gT G SL(2,R)/F : # ( gA So n H A ,c = k)) 

and 

P(K(S 0 ) = k) = /xr (gT G SL( 2 ,R)/r : # ( gA So n R A = k )). 

3.2. Expectation. To compute the expectation of the random variables K and 
EST in this setting, we use the Siegel-Veech formula 0HJ. This states that for any 
SL(2, R)-invariant measure p on R where the Siegel-Veech transform 

ns ) = E /(*) 

zeA s 

is in L 1 (p) for any / G C c (]R 2 ), there is a constant (the Siegel-Veech constant) c^ 
so that 

[ f(S) = / fdm , 

Jn J r 2 

where m is Lebesgue measure on R 2 . Applying this to our situation, we say that 
the expectations of our limiting random variables is given by a scalar multiple 
of the area of the sets H A c and R A , depending on the circle limit measure. The 
computation of Siegel-Veech constants is an active and challenging area of research, 
see, for example [T6] for seminal work. In the setting of lattice surfaces, Veech [44] 
related these constants to the covolume of SL(X, cu) in SL(2,R). 

4. Equivariant Processes 

In this section we define the axiomatic setup of equivariant measure-valued pro¬ 
cesses and state our main result Theorem im This perspective is inspired by the 
work of W. Veech [44] and J. Marklof, as well as that of A. Eskin and H. Masur Ha¬ 
lt has the great advantage that once we make the proper definitions, the proof of 
the main theorems are essentially tautologies. The power of the method lies in its 
flexibility: we will see that the axioms can be verified in several different situations. 
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4.1. Equivariant measure processes. Let n > 2, and G C GL(d,R). Let (X,g) 
denote a Borel-G-space together with a G-invariant Borel probability measure /i. 
A (G-) equivariant measure process (also known as a Siegel measure , see [33]) is a 
triple (X, g, v) where v is a map 

v : X -A A4(R d ) 

from X to the space A4(R d ) of cr-fmite Radon Borel measures on R d satisfying the 
equivariance condition 

v{gx) = g*v{x) 

for all g £ G, x £ X , where G acts linearly on R d . 

4.2. Erdos-Sziisz-Turan distributions. Given a sequence of equivariant mea¬ 
sure processes X = {(Xjv, t)n, b/v)} and a Borel subset 1Z C R d we define the 
Erdos-Sziisz-Turan distribution g = 7y(X, 1Z) on R + as the measure given by (if the 
the limit exists) 

g(X,lZ)(0,t) = lim pn(x £ A' : iz n (x)( 1Z) < t). 

N —>oo 

4.3. Equidistribution. Our main result concerns the setting where our sequence 
X = {X, g n, u}, that is, a sequence of measures g n on a fixed G-space X together 
with an assignment v. 

Theorem 4.1. Suppose g]y —> g (in the weak-* topology). Then 
g(X,lZ)(0,t) = g(x £ X : v(x)(lZ) < t). 

Proof. By assumption, our measures are all Radon Borel measures, so if t]n —> g 
in the weak-* topology, a standard approximation arguments gives that for any 
(fixed) Borel measurable subset B C X, 

Vn(B) —> g{B). 

Applying this to B = {x £ X : v{x)(lZ) < t}, we have our result. □ 

This theorem, as stated, is a tautology. The key to applying it is finding appro¬ 
priate equidistribution results that allow one to take a natural sequence X and find 
a limiting measure so that g n —> g. 

4.4. Orbits and point processes. In many of our applications, the measures gw 
will be supported on orbits of subgroups H C G and be the push-forward of some 
measure on H under the orbit map. In addition, v will often be a point process, 
that is, the assignment of a discrete set with counting measure. 

5. Equidistribution on the space of lattices 

We prove our main Diophantine results using equidistribution results for flows on 
the space of unimodular lattices. Let gu denote the Haar probability measure on 

X d = SL(d,R)/SL(d,Z). 

X d is the space of unimodular (covolume 1) lattices in R d , via the identification 

gSL{d,Z) >—> gZ d . 

Given A = gZ d £ X, we say that v £ A is primitive if 

v = gw,w £ Z d \{0}, gcd(n;) = 1 
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and denote by A pr i m the set of primitive points in A. For all of our Diophantine 
results, we will use the equivariant assignment 

gSL(d,Z)^ 6 

veqZ d ■ 

3 prim 

which, in the notation of IJD we view as a map v,i : Xj —> Af(K d ). Let m,n be 
positive integers and let d = m + n. Set 


G = SL(d,R), T = SL(d,Z), u x = ( Id Q m , H = {u x : X e Mat mx „(R)}. 

The group H is the expanding horospherical subgroup of G with respect to 

(5.1) g t = diag(e t/m ,..., e t/m , e t/n ,..., e t/n ), t > 0. 

The following Lemma is a straightforward generalisation of the the argument in the 
introduction and allows us to interpret the Erdos-Sziisz-Turan and Kesten distri¬ 
butions in terms of homogeneous dynamics. 

Lemma 5.1. Let notation be as above. Then 

EST(A, c,N) = k if and only if # (ffiog jvUxZp rim n H A , C ) = k. 

We can therefore proceed as before. Let r] N denote the measure dY on the set 
{SiogJvWi'Zprjm : 0 < ||F|| < 1}. It is well known that 

Vn — Hd, 

as N —> oo, where the convergence, as before, is in the weak-* topology. This 
seems to date to Rogers [351 P-250, (4)], who claims the result (without proof) (see 
Rogers [36j Chapter 4], for a proof of an averaged version). We refer the reader to 
Kleinbock-Margulis [2S] where a stronger statement, with a rate of convergence is 
proved. We note that, Zagier’s theorem, used in the introduction also comes with 
a rate, however the rate of convergence in these equidistribution statements does 
not shed additional light on the Erdos-Sziisz-Turan distribution. Let Xk denote the 
indicator function of the set 

{Ael : # (Aprim n Ha,c = k)} . 

The functions \k can be approximated by continuous functions with compact sup¬ 
port on the space of lattices X d+1 , so we have, as before, 


P(EST mXn (A,c,N) = k) = t] N (xk) —> V<i{xk)- 

N—yoo 


5.1. Diophantine approximation on curves. To obtain Erdos-Sziisz-Turan and 
Kesten distributions for Diophantine approximation on curves, we follow the same 
procedure above and use the following equidistribution theorem for expanding 
translates of curves due to N. Shah [ 381 . 


Theorem 5.2. Let <p : [a,b\ —> R n be an analytic curve whose image is not con¬ 
tained in a proper affine subspace. Let T be a lattice in G. Then for any Xq £ G/T 
and any bounded continuous function f on G/T, 

rb r 


(5.2) 


1 


lim 

t—► oo b — a 


f (gtu(4>(s))xo)ds = / fdg. 

Jg/t 
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6. Equidistribution ON STRATA 

6.1. Almost everywhere equidistribution. We prove Theorems 16.11 and 13.21 
using Theorem 14.11 and known equidistribution results on the space of lattices. 
Here, the equivariant assignment is given by 

S -> E A 

zEAs 

the counting measure on A§. In the notation of we denote this assignment v. 
We are interested in the Lebesgue measure of the set of 9 £ [0, 2tt ) so that 

# (rgAs (~1 Ha, c ,n) = k 

(or Ra,n )■ Applying gi og N , and using equivariance, we rewrite this as 

(A giogJvreS H Ha,c ) = k 

(respectively Ra)- 

Let t]n{S) denote the uniform probability measure M on the curve 
{ g\ og NreS : 0 <9 < 2tt} C Q g , 

and let v n = v denote the equivariant assignment. Thus, we we are in a position to 
apply Theorem l4.il with 1Z = Ha, c (or Ra)- Theorem 13. II then follows from Nevo’s 
equidistribution result which states that tin (So) -a p-u for p-u~ a.e. Sq £ R. 

6.2. Lattice surfaces. For Theorem 13.21 we restrict our universe to the subset 
SL(2,K)Sb = SL(2,R)/T, where T = SL(X 0 ,wo). Now, the sequence of measures 
can be viewed as the measures (p- supported on large circles 

{giogwreT : 0 < 9 < 27r}. 

By, for example, Dani-Smillie m , the limiting measure is the Haar probability 
measure on SL(2,R)/T, yielding our result. □ 

6.3. Other equivariant assignments. We note that there are many other equi¬ 
variant assignments (see m § 2 ] and [44] ) which can be studied in the context of 
translation surfaces. Our results, of course, apply to all such assignments. 

7. Concluding Remarks 

7.1. Complex Diophantine approximation. A complex analogue of Dirichlet’s 
theorem states that for every z G C\Q(i), there exist infinitely many p,q £ Z[z] 
such that 

\z-p/q\ < l/\q\ 2 . 

One can then consider the modified inequality | z — p/q\ < A/|g| 2 and the associ¬ 
ated Erdos-Sziisz-Turan distribution (the optimal A which allows infinitely many 
solutions for all z above is 1 /v/3, by a result of L. Ford [20]). Complex analogues of 
our results can be obtained in using analogues of our equidistribution statements. 
It would be interesting to compute the precise distribution. 

7.2. Distribution of cusp excursions. More generally, we can consider the dis¬ 
tribution of orbits non-uniform lattices A C SO(n, 1) on the Clifford plane r ra [4j. In 
this setting, the Erdos-Sziisz-Turan distribution is closely related to the distribution 
of cusp excursions for geodesic flow on HI"/A, a problem we study in |5]. 
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7.3. S-arithmetic approximation. Using our methods, one can obtain limiting 
distributions and moments in both the Erdos-Sziisz-Turan and Kesten settings for 
p-adic Diophantine approximation, and more generally in the S'-arithmetic, char¬ 
acteristic zero setting. In a forthcoming work [7], we show how to obtain such 
limiting distributions to the positive characteristic setting, proving along the way, 
equidistribution results of independent interest. 

7.4. Lines. It is natural to ask if Theorem 12.31 namely the Erdos-Sziisz-Turan and 
Kesten distributions for nondegenerate curves holds for lines. In this case, the 
distribution will depend on the Diophantine properties of the slope of the line. In 
a forthcoming work [5], we comprehensively examine this setting. 

7.5. Multiplicative approximation. The methods of this paper can be adapted 
to more general forms of the Erdos-Sziisz-Turan and Kesten distributions which 
were previously completely inaccessible. For instance, one can consider multiplica¬ 
tive and weighted forms of Diophantine approximation. In this case, one needs to 
consider equidistribution of expanding translates of more general diagonal elements. 
Let a + be the set of k tuples t = (U,..., f*,) e R fc such that 

m n 

> 0 for 1 < i < k and £* = £ tn+j ? 

*= 1 3 =1 

and for tga + define 

(7.1) g t =diag(e tl ,...,e tm ,e _tm+1 ,...,e _tfc ) 

Applying results of Kleinbock-Margulis m (see also [2Hj) which establish equidis¬ 
tribution of gt translates on SL(d + l,R)/SL(d + 1,Z), it should be possible to 
obtain distributional results in this setting. Similarly, using a result of Shah 551. 
one can also obtain Erdos-Sziisz-Turan and Kesten distributions in the multiplica¬ 
tive setting for curves. 

7.6. Inhomogeneous approximation and affine lattices. Our methods can 
also be applied to the set of affine lattices AXd = SL(d, R) k R d / SL(d, Z) x Z d to 
study the distribution of points in affine lattices, with applications to inhomoge¬ 
neous approximation, see m for related results. 

7.7. Cut-and-project quasicrystals. The fine scale statistics of cut-and-project 
quasicrystals have been studied numerically by Baake-Gotze-Huck-Jacobi [5] and 
using methods of homogeneous dynamics by Marklof-Strombergsson [32]. Following 
the approach of [52] . our methods should extend to this setting, and it would be 
very interesting to work out the explicit distributions. 
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